On the offset of the DM cusp and the interpretation of the 130 GeV line as a DM 

signal 



Dmitry Gorbunov 

Institute for Nuclear Research of Russian Academy of Sciences, 
60-th October Anniversary pr. 7a, 117312 Moscow, Russia and 
Moscow Institute of Physics and Technology, 
Institutsky per. 9, 141 700 Dolgoprudny, Russia 



(N 

o 

O 

Q 



O 

U 



> 
oo 

00 

o 



X 



Peter Tinyakov 

Service de Physique Theorique, Universite Libre de Bruxelles (ULB), 
CP225 Boulevard du Triomphe, B-1050 Bruxelles, Belgium 

We show that the cusp in the dark matter (DM) distribution required to explain the recently 
found excess in the gamma-ray spectrum at energies ~ 130 GeV in terms of the DM annihilations 
cannot survive the tidal forces if it is offset by ~ 1.5° from the Galactic center as suggested by 
observations. 



I. INTRODUCTION 



Recently, a line-like feature in the gamma-ray spec- 
trum observed by the Fermi LAT experiment in the di- 
rection of the Galactic center has been found. It has been 
suggested that it can be interpreted in terms of the dark 
matter (DM) annihilation signal [l|, 0] ■ 

More refined analysis has confirmed the significance 
of the feature in the spectrum d, 0]. The best fit to 
the DM annihilation line was obtained with the modified 
Navarro-Frenk- White (NFW) DM profile with the inner 
slope of a = —1.2. It also revealed the offset of the sig- 
nal with respect to the Galactic center by about 1.5°, 
which corresponds to the projected distance of ~ 200 pc. 
The existence of the offset, if confirmed, raises a ques- 
tion about the viability of the DM interpretation of the 
observed gamma-ray excess. 

In order to address this question, numerical simula- 
tions of the Galactic center including the effect of the bar 
were used in Ref. [e$ . It was found that an offset of a few 
hundred parsec between the GC and the maximum of the 
DM distribution could in principle exist (see also Ref. [6] 
and references therein). However, the DM distribution in 
this case is not cuspy but cored and has an overdensity 
with respect to the central region of only about 10 — 20%. 
In addition, it was found that the core density is formed 
by the gravitationally unbound DM particles. 

Here we address the same question from a different 
perspective. Making use of the analytical estimates we 
show that only a small central region of the DM cusp can 
survive the tidal forces produced by the baryons which 
dominate the gravitational potential near the Galactic 
center. The surviving part of the cusp is too small to 
explain the observed feature in the 7-ray spectrum. 



II. BARYONIC AND DM PROFILES IN THE 
GALACTIC CENTER 

The baryon distribution in the GC is known from the 
2 fim light distribution (see, e.g., Refs. [7J, l8|) and con- 
firmed by the study of the kinematic properties of the 
OH/IR stars @. 

For our purposes the exact behavior of the baryon den- 
sity is not necessary, and a crude approximation in the 
inner ~ 200 pc is sufficient. From Fig. 10 of Ref. @ we 
adopt the following approximation for the baryonic mass 
Afg(r) enclosed within the radius r, 



M B (r) = 6 x 1O 8 M 



100 pc 



1.25 



(1) 



where Mq is the solar mass. From this equation, the 
baryonic density pb(t) is 



PB (r) = 60M Q pc- 



100 pc 
= 2.3 x 10 3 GeVcm" 3 



-1.75 



100 pc 



-1.75 



This relation is consistent, within the errors, with the 
one derived in Refs. It is clear from eq. (P) that 

at distances > 100 pc from the Galactic center the effect 
of the central black hole is subdominant, and we ignore 
it in what follows. 

Now we turn to the DM distribution. In Ref. [4f sev- 
eral DM profiles were shown to fit the gamma-ray data, 
namely the Einasto profile [ll[ and modified NFW pro- 
files of the form 



p(r) 



(r/r s ) a {l+r/r s ) 3 



(2) 



where r s — 20 kpc and the slope a ranges form a = 1 to 
a = 1.3, with the best fit value a = 1.2. The normaliza- 
tion factor p s — 0.27 GeV/cm 3 is fixed by setting the DM 
density around the Earth to 0.4GeV/cm 3 . To illustrate 
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our point it is sufficient to consider the simpler and more 
cuspy NFW profile ((2J. The less cuspy Einasto profile is 
subject to even stronger tidal effects. 

As a first step, let us find the size of the region which 
is responsible for the gamma-ray signal assuming the lat- 
ter is produced by the DM annihilations. The observed 
signal corresponds to the luminosity [J] 

L Q = (1.7 ± 0.4) x 10 36 photons/s. 

On the other hand, the luminosity of a spherical region of 
the size r centered on the DM distribution can be written 
as follows, 



L(r) 



'DM 



[av) p 2 s I{r/r s ) 



(3) 



where (av) ~ 2 x 10 _27 cm 3 /s is the velocity-averaged 
DM annihilation cross section Q , Wdm — 140 GeV is 
the DM mass and 



I(r/r s ) = 



r jr s 



».2u 



(1+x) 



6-2a 



(4) 



Equating this to the observed luminosity gives the fol- 
lowing equation for the size r of the emission region, 



I(r/r s ) 



DM 



47T (av) p 2 s rf 



= 7.7 x 10" 



This equation is easily solved by noting that the solution 
corresponds to small values of r/r s for which the integral 
in eq. Q can be simplified, 



I(r/r s ) 



■{r/rs) 



3~2a 



3 - 2a 

For the best fit case a — 1.2 this gives 
r ~ 0.006 r s = 120 pc. 



(5) 



Note that when obtaining this estimate we have assumed 
(cf. eq. ©) that all of the DM mass is converted into 
photons, as in the case of the annihilation into 77. If the 
efficiency were lower, as it would be in the case of the 
annihilation into a single- photon final state P, the 
size of the contributing region would be even larger. 

It is instructive to calculate the total amount of DM 
contained within a given distance from the cusp. Making 
use of the relation r/r s <C 1 one finds 



1 



(r/ r s 



MdmW = 471-^ r 3 



where r is the distance form the center of the cusp. This 
implies for a = 1.2 



M DM (r) = 2.9 x 10 7 M G 



100 pc 



(6) 



The latter value has to be compared to eq. (JTJ . Clearly, at 
distances ~ 200 pc from the Galactic center the baryons 
give a dominant contribution to the total mass. Thus, 
they dominate the gravitational potential except in the 
vicinity of the cusp. 



III. TIDAL STRIPPING OF THE OFFSET DM 
CUSP 

The DM cusp that is offset with respect to the bary- 
onic distribution is subject to tidal forces. To estimate 
the importance of these forces one may compare the dif- 
ference of the gravitational pull of baryons at different 
parts of the DM distribution and the gravitational force 
from the DM itself. For the cusp to survive the gravi- 
tational force from baryons has to be smaller than the 
force from DM. This leads to the condition 

GM DM (rf) GM B {r) J 



(P 



(7) 



where r ~ 200 pc is the offset distance and it was assumed 
that d <C r. Making use of eqs. (JTJ and (JBJ one finds for 
a = 1.2, 



d < 22 pc 



(8) 



Thus, the tidal radius (the maximum radius where DM 
particles survive the stripping by tidal forces) is much 
smaller than the offset distance, which justifies the ap- 
proximation used. 

More accurately, the tidal radius d can be calculated 
by making use of the formalism developed in Ref. [l2j | . In 
principle, one should distinguish tidal radii corresponding 
to prograde, radial and retrograde DM orbits. However, 
in the long time limit (that is, at times much longer than 
the period of the orbital motion ~ 10 7 yr) these converge 
to the smallest of the three. Making use of eq. (21) of 
Ref. 113. one finds at a = 1.2 



5pc 



(9) 



At smaller values of a the DM cusp is weaker and the 
tidal radius is slightly smaller, while at larger values of 
a it is slightly larger, always being of the same order as 
given by eq. ([S]). Thus, only a very small central part of 
the cusp can survive the tidal disruption. This part is 
much smaller than the region responsible for the annihi- 
lation signal. 



IV. CONCLUSIONS 

From the mismatch between eqs. ([5]) and ([9]) it is clear 
that the part of the cusp that can survive the tidal strip- 
ping is insufficient to explain the DM signal. First, its 
angular size is about 2'. Such a source has an angular size 
smaller than the Fermi-LAT point-spread function [l3| 
and would appear as point-like, which is not compatible 
with the morphology of the observed excess. 

Second, according to eq. ([3]), the annihilation signal 
from such a small region would be reduced by a factor 
of about ~ 7 if the annihilation cross section of (av) ~ 
2 x 10 -27 cm 3 /s is assumed. To make the signal from the 
surviving part of the cusp compatible with observations, 
one would have to increase the cross section by the same 
factor, which would be in contradiction with the limits 
from Fermi LAT [H. 
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